Coiling and Supercoiling of Vortex Filaments in Oscillatory Media 
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We study the behavior of vortex filaments subject to a uniform density of phase twist in oscillatory 
media described by the complex Ginzburg-Landau equation. The first instability is a supercritical 
Ifopf bifurcation to stable propagating helical vortices. The secondary instability, also a forward 
Ifopf bifurcation, leads to quasiperiodic supercoiled filaments. The structural changes undergone 
by these dynamical objects are akin to those of twisted elastic rods, in spite of the presence of the 
ribbon component of the twist, particular to phase singularities. 



Spontaneous oscillations in nonlinear extended media 
are commonly observed in nature. Examples include 
reaction-diffusion systems certain regimes of fluid 
flows 1^], and biological systems with intrinsic clocks 
The complex Ginzburg-Landau equation (CGLE) 



dtA = A+{l + ia)V^A - (1 + i(3)\A\'^A 



(1) 



where ^ is a complex amplitude describing the slow mod- 
ulations of the oscillations and a and (3 are two real pa- 
rameters, plays a central role in the analysis of the be- 
havior of oscillatory media, especially since its validity 
extends far beyond the vicinity of the Hopf bifurcation 
where it can be formally derived, 

The CGLE exhibits the localized structures that of- 
ten determine the spatiotemporal organisation of the 
medium. In two dimensions, A — Q generically at point 
vortices which may be sources of spiral waves. In three 
space dimensions, the topological defects are filaments 
around which "scroll waves" can be emitted. The de- 
fects are phase singularities, determined by the topology 
of the phase field aigA. Consequently, they are not de- 
scribed solely by their geometrical structure as given by 
the Frenet frame at each point along the filament; one 
must also specify how the phase field itself is twisted 
along the vortex (ribbon twist). ||^ 

As in lower space dimensions, the structure and dy- 
namics of three-dimensional (3D) phase singularities is 
best studied in the parameter region where the scroll 
waves are stable. Recent simulations and analytical stud- 
ies in this region have shown that initially prepared 
untwisted vortex rings shrink and disappear. Disordered 
states [D and helical structures arising from sponta- 
neous stretching and bending of untwisted vortex fila- 
ments have been found in this region. 

In this Letter, we investigate numerically the effects of 
a finite uniform density 7 of ribbon twist per unit length 
on the simplest configuration, a straight infinite filament 
in the regime where both the emitted scroll waves and 
the 2D spiral cores are stable. While the limit of small 
twist was discussed in [Q, a detailed study of the effects 
of twist has not been carried out for the CGLE, even 



though twisted filaments are easily found in the chaotic 
regimes of the 3D CGLE and have also been observed 
in experiments on the Belousov-Zhabotinsky (BZ) reac- 
tion in the excitable regime. We find that there is a 
finite twist density 7c (a, (3) beyond which the infinite fil- 
ament undergoes a supercritical Hopf bifurcation which 
saturates to produce helical vortices. We also charac- 
terize the secondary instability, which takes the form of 
another Hopf bifurcation leading to the supercoiling of 
the primary helices. These two bifurcations are analyzed 
from a topological viewpoint. In particular, the writhe 
is a linear function of the excess twist after the first in- 
stability, and the ribbon twist also varies linearly. This 
suggests analogies between the dynamical structures ob- 
served and the behavior of twisted elastic rods. 

Our results extend similar early observations in ex- 
citable media. Of special relevance to the present 
work are studies of helical solutions of various reaction- 
diffusion models which have shown that an ini- 
tially straight vortex filament may adopt a helical form 
with fixed radius for large enough twist. Helical filaments 
have also been observed in 3D magnetic resonance images 
of the BZ reaction. We believe our findings apply to 
the excitable case and put those results in a more general 
setting by elucidating the dynamical nature of the insta- 
bilities. Our work should also lead to an understanding 
of the behavior of more general twisted vortex lines. 

The spiral solution of the 2D CGLE is important for 
the study of 3D vortex lines: cross sections of scroll waves 
are spirals, up to correction terms calculated in 0. The 
spiral solution is given in polar coordinates (r, 9) by 



A^{r, 6, t) = F{r) exp i [t^^t ±61 + V(r)] 



(2) 



with ijjs = — /3 + fcg (/3 — a), where F{r) and ip{r) are two 
real functions whose asymptotic behavior for r — > and 
r — > 00 is known, and kg is the asymptotic wave number. 
p6[ In spite of the absence of a general explicit form, 
many properties of this solution are known, in particular 
in the parameter region where the asymptotic emitted 
planewave and the discrete core modes are linearly sta- 
ble. We first consider straight twisted filaments in 
this domain of the (a, /3) plane and use the 2D stable 
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spiral solution to generate initial conditions efficiently by 
shifting the phase of the spiral regularly along the third 
spatial dimension z to create twist: 9 6 + ^z. In the 
numerical simulations reported below, such initial condi- 
tions have been used in finite boxes of "height" Lz with 
periodic boundary conditions in z, so that the total twist 
is a multiple of 27r: 7L2 — 27rnt. 

The link number Lk, total twist Tw and writhe Wr of 
the filament ribbon obey the following law p^,p[|: 



Lk = Tw + Wr = Twf + Tw,. + Wr 



(3) 



where Twf — § T{s)ds and Tw, — ^(dip/ds) ds are, re- 
spectively, the total Frenet and ribbon components of 
the twist, defined here in terms of local quantities cal- 
culated in the Frenet frame (t, n, b) parametrized by the 
arclength s, the torsion r and the angle (p between a given 
ribbon curve and the binormal vector b. For closed 3D 
ribbons, the link number Lk is invariant j^]. In finite 
boxes periodic in z, it is conserved if the ribbon is not 
"cut" , i.e. no reconnection occurs, either of the filament 
with itself or with another filament, and we have, simply: 



Lk = 2tt nt = 



(4) 



The constraint imposed by the conservation law has im- 
portant consequences since it implies that filaments with 
different 7 values are not continuously related; thus, one 
must consider 7 as an independent parameter. 
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1. Stable helical vortex for a = 1.7, /3 = 0, 7 = 0.098 
= 512 {Lx = Ly = 128). Top: perspective view of the 
3D isosurface \A\ = 0.6 colored by the phase field (27r-periodic 
grey scale). Note the presence of ribbon twist (7h = 0.196, 
Twr/Lz = Bottom left: {x,z) projection of the filament 

at t = 2495 (dashed line) and t = 2500 (solid line). Bottom 
right: growth of {R)z (solid line) and time series of the x co- 
ordinate (dashed line) of a point on the filament at a given z. 

Consider a single twisted filament, generated as de- 
scribed above, evolving under the CGLE. Our simula- 
tions employed no-flux boundary conditions in the x and 
y directions while the box was periodic in z, the axis of 
the filament. All numerical results have been checked by 
changing the box size and the spatiotemporal mesh size. 

For small 7, the solution quickly relaxes to take into 
account the wavenumber and frequency shifts between 



the 2D spiral waves and the 3D scroll waves. This (sta- 
ble) solution of the 3D CGLE is Af (r, 6*, z, t): 

Af = Vl " l'^F{r') exp ziwft ±{9 + jz) + ij{r')] , (5) 



with 



where r' = r-\/l — 7^ and Uf = u!s{i — 7^) 
F, and given by the 2D spiral solution (y). We have 
confirmed that our numerical solutions coincide with the 
above exact solution when there is no core instability. 
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FIG. 2. a = 1.7 and P = Q: R{j) (left) and Wr(7)/L^ 
(right) (7c ~ 0.032). Insets: R^ij) (left) and 7h(7) (right). 
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FIG. 3. Left: phase diagram of the 3D CGLE. The re- 
gion of interest lies between the 2D spiral core instability 
(dash-dotted) line, the phase instability l+af^ = (solid line), 
and the a — P line (dashed) where the spiral core is infinite. 
The thin lines are the 7c = 0.049 (circles) and 7c = (squares, 
from [^) level curves of the 7c(q,/3) surface. Right: 7c(a) for 
/3 = 0. A quadratic fit gives 7c = for a = 1.95(15). 

For large enough 7 this relaxation is followed by the 
growth of periodic modulations along the filament, which 
becomes a helix of pitch Ah and well-defined radius R 
(Fig. p. The helix rotates with constant angular veloc- 
ity Lu and any point on the helix moves along its axis at 
a constant velocity v = uj/^h with 7h = 27r/Ah. These 
results indicate that the straight filament has undergone 
a Hopf bifurcation to traveling waves. The direction of 
propagation and rotation is determined by the sign of 
7, which breaks the z ^ —z parity symmetry. Varying 
7 continuously across the instability threshold 7c, the 
behavior expected for a supercritical bifurcation is ob- 
served: ~ (7 ~ 7c), while v and Ah vary smoothly 
(Fig. The forward nature of the bifurcation is con- 
firmed by the growth transient (Fig. |l| bottom right). 
Similarly, keeping 7 constant and varying (a, f3) toward 
the 2D core-instability line, a bifurcation point is passed 
and the expected small- i? behavior (e.g. R^ ~ (a — ac) at 
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fixed /?) is observed (not shown). This suggests the exis- 
tence of a critical surface 7c (a, P) deUmiting the stabihty 
of straight twisted filaments. Our simulations yielded the 
approximate location of the 7c = 0.049 line in the (a,/9) 
plane (Fig. ||lcft). The complete 7c(ck,/3) surface could 
be determined from a linear stability analysis of solution 
(^, possibly along the lines of a rather difficult task 
left for future work [|o|. 

The 7c = line determines the stability domain of the 
untwisted filament. In it was argued that this solu- 
tion, which is just the trivial extension of the 2D spiral, 
is unstable in the region of the 2D core instability. It 
was further argued recently that the instability region of 
the untwisted filament extends beyond the 2D core line 
. To attack this problem from the viewpoint of twisted 
filaments, we have performed simulations for /? = to de- 
termine 7c (a) (Fig. 1^ right). We find, by extrapolation, 
7c = for a = 1.95(15), in agreement with Our re- 
sults indicate that the instability of the 7 = filaments 
does not lead to helices. In this special case, the z — > —z 
parity symmetry is preserved and one expects, a priori, 
the superposition of opposite traveling waves at the linear 
stage. At the nonlinear level, we observe the saturation 
of the two corresponding amplitudes at the same value 
and the appearance of a flat, z-periodic filament rotating 
uniformly around the z axis, the equivalent of standing 
waves in the amplitude representation (Fig. Q) . 




FIG. 4. Asymptotic state of an unstable untwisted filament 
{a — 2.5, /3=0). Top: 3D isosurface =0.6 colored by the 
phase field. Bottom: space-time plot {t running from left 
to right) of the phase ^(z,*) = a,Tcta.n{y{{z) /x{{z)) where 
{x{{z),y{{z), z) parametrize the filament. (27r-periodic grey 
scales 0=black, 7r=white, were used). 

We now consider our results at non-zero twist from 
a topological viewpoint. Since no reconncction occurs 
across the transition, Lk remains constant. Given Ah 
and R, all terms of conservation law (||) can be calcu- 
lated easily. For the straight filament, the Frenet frame 
is degenerate; Twf and Twj. cannot be distinguished. For 
the stable helix, on the other hand, the ribbon twist is 
given by Twi/Lk = 1 — 7h/7. Moreover, k and r are 
constant along the z axis: K/{R'yh) = t = + R'^j^). 

We then have [||: 



Twf _ 7h _ Wr _ 7h 

Lk - 7 Lk " j^TTW^ ■ ^ ' 

Exploring a large part of the (a, /?) domain of interest, 
we consistently found that Ah is larger than As, the wave- 
length of the 2D spiral, but of the same order, and that 
the two wavelengths vary similarly with (a, (3) [l0| . Thus 
the typical scale of the 3D instability depends essentially 
on the wavenumber selected by 2D spiral solution. We 
also find that 7h > 7, which implies that Twj. < 0. The 
ribbon twist "compensates" the strong torsion due to the 
instability. For {a = 1-7,(3 = 0), Ah decreases with in- 
creasing 7 > 7c (Fig. |[ right), and thus the magnitude 
of the Twr increases. In fact, our data indicates a lin- 
ear variation with 7, but we cannot rule out a weakly 
quadratic dependence. Moreover, in agreement with (|^) 
for small i?^ ~ 7 — 7c, Wr ~ (7 — 7c), but, remarkably, 
this behavior also extends far from 7c (Fig. ||). (Again, a 
weakly quadratic dependence is still possible.) Thus, all 
topological quantities seem to vary linearly above thresh- 
old, and, in particular, a constant proportion of the excess 
twist 7 — 7c is converted into writhe. This is akin to the 
behavior of elastic rods |||,|2^ . 

Moving away from the 7c(a,/3) surface, the filament 
undergoes a secondary instability at some critical value 
7^ > 7c: a long- wavelength modulation appears on the 
stable regular helix, eventually producing a supercoiled 
filament (Fig. ^). No reconncction event is observed. 
This structure is characterized by two well-defined wave- 
lengths, Ai and A2, with Ai ~ Ah and A2 3> Ah (Fig. ||, 
middle left). In general, for an infinite filament, the two 
wavelengths are incommensurate and the superhelix is a 
quasiperiodic object whose projection on the {x, y) plane 
produces a characteristic "flower" pattern (middle right). 
Associated with these wavelengths are two frequencies, or 
two velocities vi and v-2 (middle left). 

The secondary instability is also a supercritical Hopf 
bifurcation. Its forward nature and the nonlinear sat- 
uration are apparent in the time evolution of an ini- 
tially straight filament, which quickly turns into a weakly 
unstable helix before bifurcating toward a superhelix 
(Fig. ^, bottom left). Varying 7, the mean radius (R)z 
is continuous but not differentiable at 7^, and the am- 
plitude of the secondary modulations can be measured 
by (R)^ — Rh, where i?h is the radius of the (unstable) 
primary helix (measured, e.g., from the transient helix 
stage) (Fig. || bottom right). As expected for a forward 
Hopf bifurcation, near threshold our data is consistent 
with {{R)z - Rh)^ (7 - 7c)- Even though the ampli- 
tude of this secondary instability saturates, we cannot 
rule out a chaotic asymptotic behavior. 

Our measurements (from Fourier spectra) show that 
1/Ah = 1/Ai + I/A2, at least for 7 — 7c not too large. 
Since the link number Lk is constant (no reconncction 
occurred), this relation implies that the ribbon compo- 
nent of the twist Twr varies smoothly across the tran- 
sition. This is corroborated by direct measurements 
of Twr and the other topological quantities (Wr, Twf) 
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across the secondary instability. During the evolution 
shown in Fig. I (bottom left), Wi{t) is constant through 
the second morphological change (bottom left) since the 
mean torsion remains constant; r(s), however, becomes 
Ai-periodic. Finally, the available data does not allow to 
decide whether these quantities still vary linearly with 7. 




phase mode, it can be seen as accomodating the torsion 
and writhe of the bifurcated structure. Indeed, although 
it is a dynamical object supporting propagating waves, 
our filament coils and supercoils under twist like an elas- 
tic rod or a DNA molecule (see [| 20 and references 
therein). Moreover, the simple behavior of Wr supports 
the adoption a topological viewpoint when studying the 
dynamics of 3D phase singularities. We note, however, 
that the ribbon twist might play a more important role 
in excitable media or systems with complex oscillations, 
for which the gauge invariance of the phase is lost. 
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FIG. 5. Secondary Hopf bifurcation for a = 1 and 
P = —0.5. This data was obtained under the contraint that 
Twr = 0, a small perturbation here, since 7^ ~ 7 for these 
parameter values. Top: 7 = 0.27; 3D isosurface \A\ = 0.6 
colored by the phase field (27r-periodic grey scale); Middle 
left: space-time plot of 4>{z,t) (see Fig. 0). Portion of size 
L = 600 of a box of length = 2048 shown during At = 1825 
with an initial straight filament [Lx — Ly = 92, 7 — 0.295). 
27r-periodic grey scale (0 =black, n =white), except for the 
points where dz(f) < which are plotted in black regardless 
of the value of (f). This enhances the "petals" , or small loops, 
in the superhelix, since only points of these loops have neg- 
ative dz(l}- For the primary helix, observed as a transient, 
dz(j) = 'yh > everywhere. Wavelengths Ah, Ai, and A2 and 
velocities v, vi and V2 are indicated. Middle right: projection 
of superhelix on the {x,y) plane. The fiower pattern is peri- 
odic due to the finite box length. Bottom left: 7 = 0.27: time 
series of (R)^ and an its rms over z (left scale), Wr/L^ and 
Gt the rms of the torsion (right scale). Bottom right: (Ii}z{"/) 
across the two bifurcation points 7c — 0.171 and 7^ — 0.267. 
Dashed line: _Rh for the (weakly) unstable helices. 

Although the instabilities studied here are rather tra- 
ditional Hopf bifurcations, they are topologically con- 
strained. The existence of the ribbon component of the 
twist, particular to the filaments studied here, comes into 
play in a somewhat "passive" way. Associated to the soft 



[1] Chemical Waves and Patterns, eds. R. Kapral and K. 

Showalter, (Kluwer, Dordrecht, 1994). 
[2] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 

851 (1993). 

[3] A. Goldbeter, Biochemical Oscillations and Cellular 
Rhythms, (Cambridge U. P., Cambridge, 1996); A. T. 
Winfree, When Time Breaks Down, (Princeton U. P., 
Princeton, NJ, 1987). 

[4] Y. Kuramoto, Chemical Oscillations, Waves and Turbu- 
lence, (Springer, Berlin, 1984); A. C. Newell, Envelope 
Equations, (Am. Math. Soc, Providence, RI, 1974). 

[5] For recent experimental examples, see, e.g., T. Leweke 
and M. Provansal, Phys. Rev. Lett. 72, 3174 (1994); Q. 
Ouyang and J.-M. Flesselles, Nature 379, 6561 (1996). 

[6] M. Tabor and I. Klapper, Nonlin. Sc. Today, 4:1, 7 
(1994); 4:2, 12 (1994). 

[7] M. Gabbay, E. Ott and P. N. Guzdar, Phys. Rev. Lett. 
78, 2012 (1997). 

[8] I. S. Aranson and A. R. Bishop, Phys. Rev. Lett. 79, 
4174 (1997). 

[9] I. S. Aranson, A . R. Bishop and L . Kramer, xxx.lanl.gov 

e-Print archive, t3att-sol/9712004 , 
[10] G. Rousseau, H. Chate, and R. Kapral, in preparation. 
[11] A. M. Pertsov, R. R. Aliev and V. I. Krinsky, Nature 

345, 419 (1990). 
[12] For a review see, A. T. Winfree, SIAM Rev. 32, 1 (1990); 

in §, p. 1. 

[13] C. Henze, E. Lugosi and A. T. Winfree, Can. J. Phys. 
68, 683 (1990). 

[14] J. P. Keener, Physica D 269 (1988); J. J. Tyson and J. 
P. Keener, in [Q], p. 93. 

[15] A. Cross, The Development of 3D MRI of the BZ Reac- 
tion, Ph.D. thesis, U. of New Brunswick, 1998. 

[16] P.S. Hagan, SIAM J. Appl. Math. 42, 762 (1982). 

[17] I.S. Aranson, L. Kramer, and A. Weber, Phys. Rev. Lett. 
72, 2316 (1994). 

[18] J. H. White, Am. J. Math. 91, 693 (1969). 

[19] With no- flux conditions in z, the link would not be a 
conserved quantity, complicating further the picture. 

[20] A. Goriely and M. Tabor, Phys. Rev. Lett. 77, 3537 
(1996); Physica D 105, 20 (1997); ibid. 105, 45 (1997); 
preprint (1997), to appear in Proc. Roy. Soc. A. 



4 



